At its nearest, northern Australia is just over 400 km from an active convergent plate margin. This complex and unique tectonic region combines active subduction and the collision of the Sunda-Banda Arc with the Precambrian North Australian Craton (NAC) near the Timor Trough and continues through to the New Guinea Highlands. Ground-motions generated from earthquakes on these structures have particular significance for northern Australian communities and infrastructure projects, with several large earthquakes in the Banda Arc region having caused ground-shaking-related damage in the northern Australian city of Darwin over the historical period. There are very few, if any, present-day tectonic analogs where cold cratonic crust abuts a convergent tectonic margin with subduction and continent-continent collision. Ground motions recorded from earthquakes in typical subduction environments are highly attenuated as they travel through young sediments associated with forearc accretionary prisms and volcanic back-arc regions. In contrast, seismic energy from earthquakes in the northern Australian plate margin region are efficiently channelled through the low-attenuation NAC, which acts as a waveguide for high-frequency earthquake shaking. As such, it is difficult to select models appropriate to the region for seismic hazard assessments. The development of a far-field ground-motion model to support future seismic hazard assessments for northern Australia is discussed. In general, the new model predicts larger ground motions in Australia from plate margin sources than models used for the 2018 National Seismic Hazard Assessment of Australia, none of which were considered fully appropriate for the tectonic environment. Short-period ground motions are strongly dependent on hypocentral depth and are significantly higher than predictions from commonly-used intraslab ground-motion models at comparable distances. The depth dependence in ground motion diminishes with increasing spectra periods. <b>Cite this article as</b> Allen, T. I. (2021). A Far-Field Ground-Motion Model for the North Australian Craton from Plate-Margin Earthquakes, <i>Bull. Seismol. Soc. Am. </i><b> 112</b>, 1041–1059, doi: 10.1785/0120210191

Seismic hazard models, commonly produced through probabilistic seismic hazard analysis, are used to establish earthquake loading requirements for the built environment. However, there is considerable uncertainty in developing seismic hazard models, which require assumptions on seismicity rates and ground-motion models (GMMs) based on the best evidence available to hazard analysts. This paper explores several area-based tests of long-term seismic hazard forecasts for the Australian continent. ShakeMaps are calculated for all earthquakes of MW 4.25 and greater within approximately 200 km of the Australian coastline using the observed seismicity in the past 50 years (1970-2019). A “composite ShakeMap” is generated that extracts the maximum peak ground acceleration “observed” in this 50-year period for any site within the continent. The fractional exceedance area of this composite map is compared with four generations of Australian seismic hazard maps for a 10% probability of exceedance in 50 years (~1/500 annual exceedance probability) developed since 1990. In general, all these seismic hazard models appear to be conservative relative to the observed ground motions that are estimated to have occurred in the last 50 years. To explore aspects of possible prejudice in this study, the variability in ground-motion exceedance was explored using the Next Generation Attenuation-East GMMs developed for the central and eastern United States. The sensitivity of these results is also tested with the interjection of a rare scenario earthquake with an expected regional recurrence of approximately 5,000 - 10,000 years. While these analyses do not provide a robust assessment of the performance of the candidate seismic hazard for any given location, they do provide—to the first order—a guide to the performance of the respective maps at a continental scale. This paper was presented at the Australian Earthquake Engineering Society 2021 Virtual Conference, Nov 25 – 26.

A database of recordings from moderate-to-large magnitude earthquakes is compiled for earthquakes in western and central Australia. Data are mainly recorded by Australian National Seismograph Network (ANSN), complemented with data from temporary deployments, and covering the period of 1990 to 2019. The dataset currently contains 1497 earthquake recordings from 164 earthquakes with magnitudes from MW 2.5 to 6.1, and hypocentral distances up to 1500 km. The time-series data are consistently processed to correct for the instrument response and to reduce the effect of background noise. A range of ground-motion parameters in the time and frequency domains are calculated and stored in the database. Numerous near-source recordings exceed peak accelerations of 0.10 g and range up to 0.66 g, while the maximum peak velocity of the dataset exceeds 27 cm/s. In addition to its utility for engineering design, the dataset compiled herein will improve characterisation of ground-motion attenuation in the region and will provide an excellent supplement to ground-motion datasets collected in analogue seismotectonic regions worldwide. This paper was presented at the Australian Earthquake Engineering Society 2021 Virtual Conference, Nov 25 – 26.

Public concerns have been raised about the potential for induced seismicity as state and territory governments lift moratoriums on hydraulic stimulation activities for the exploration and extraction of unconventional hydrocarbons. The Scientific Inquiry into Hydraulic Fracturing in the Northern Territory articulated the need for a traffic-light system “to minimise the risk of occurrence of seismic events during hydraulic fracturing operations” within the Beetaloo Sub-basin. A temporary seismic network (Phase 1) was deployed in late 2019 to monitor baseline seismic activity in the basin. Based on the data analysed herein (November 2019 – April 2021), no seismic events were identified within the area of interest suggesting that the Beetaloo Sub-basin is largely aseismic. Observations to date indicate that there is potential to identify events smaller than ML=1.5 within the basin. The recent installation of ten semi-permanent stations for continuous real-time monitoring will contribute to ongoing baseline monitoring efforts and support the implementation of an induced seismicity traffic-light system. The outcome of this study will be used to build knowledge about potential human-induced seismic activity in the region that may be associated with unconventional hydrocarbon recovery. This paper was presented at the Australian Earthquake Engineering Society 2021 Virtual Conference, Nov 25 – 26.

Here we undertake a statistical analysis of local magnitudes (ML) calculated using the two real-time earthquake monitoring software platforms use by Geoscience Australia (GA) since 2005, Antelope and Seiscomp. We examine a database of just over 10 years duration, during a period in which both systems were in operation and over 4000 earthquakes were located and magnitudes estimated. We examine the consistency of both single-station and network ML estimates of both systems, with a view toward determining guidelines for combining them into a single catalogue, as well as for determining best practice in the for the estimation of local magnitudes for regions of sparse seismic networks. Once this guidance has been developed, it is the intention of GA to re-process magnitudes for all earthquakes using a consistent approach where digital data are available and can be integrated within the currently-used SeisComP system. This paper was presented at the Australian Earthquake Engineering Society 2021 Virtual Conference, Nov 25 – 26.

The Mwp 6.1 Petermann Ranges earthquake that occurred on 20 May, 2016 in the Central Ranges, NT, is the largest onshore earthquake to be recorded in Australia since the 1988 Tennant Creek sequence. While geodetic and geophysical analyses have characterized the extent of surface rupture and faulting mechanism respectively, a comprehensive aftershock characterization has yet to be performed. Data has been acquired from a 12-station temporary seismic network deployed jointly by the ANU and Geoscience Australia (GA), collected from five days following the mainshock to early October. Taking advantage of enhanced automatic detection techniques using the SeisComP3 real-time earthquake monitoring software within the National Earthquake Alerts Centre (NEAC) at GA, we have developed a comprehensive earthquake catalogue for this mainshock-aftershock sequence. Utilising the NonLinLoc location algorithm combined with a Tennant Creek-derived velocity model, we have preliminarily located over 5,800 aftershocks. With additional spatio-temporal analyses and event relocation, our objective will be to use these aftershocks to help delineate the geometry of the headwall rupture along the Woodroffe Thrust. These high-resolution aftershock detection techniques are intended to be implemented in real-time within the NEAC following future significant Australian intraplate earthquakes. This paper was presented at the Australian Earthquake Engineering Society 2021 Virtual Conference, Nov 25 – 26.

The 22 September 2021 (AEST) moment magnitude MW 5.9 Woods Point earthquake was the largest in the state of Victoria’s recorded history. The ground motions were felt throughout the state of Victoria and into neighbouring states New South Wales and South Australia. Minor damage was reported in the city of Melbourne and in some regional towns close to the epicentre. This event was captured on many high-quality recorders from multiple sources, including private, university, and public stations. These recordings provide a rare opportunity to test the validity of some ground motion models thought to be applicable to the southeast region of Australia. This paper presents spectral acceleration and attenuation comparisons of the Woods Point earthquake event to some ground motion models. The results of this paper provide further evidence that the attenuation characteristics of southeastern Australia may be similar to that in central and eastern United States, particularly at shorter distances to the epicentre. This paper was presented at the Australian Earthquake Engineering Society 2021 Virtual Conference, Nov 25 – 26.

Seismic hazard modelling is a multi-disciplinary science that aims to forecast earthquake occurrence and its resultant ground shaking. Such seismic hazard models consist of a probabilistic framework that models the flow of uncertainty across a complex system; typically, this includes at least two model-components developed from earth science: seismic source models, and ground motion prediction models. Although there is no scientific prescription for the length of the forecasting time-window, the most common probabilistic seismic hazard analyses (PSHA hereafter) consider forecasting probabilities of ground shaking in time windows of 30 to 50 years. These types of models are the target of this review paper. Although the core methods and assumptions of such a modelling have largely remained unchanged since they were first developed more than 50 years ago, we will review the most recent initiatives which are facing the difficult task of meeting both the increasingly sophisticated demands of society and keeping pace with advances in our scientific understanding. A need for more accurate and precise hazard forecasting must be balanced with increased quantification of uncertainty and new challenges such as moving from time-independent hazard to forecasts that are time-dependent and specific to the time-period of interest. Meeting these challenges requires the development of science-driven models which integrate at best all information available, the adoption of proper mathematical frameworks to quantify the different types of uncertainties in the source and ground motion components of the hazard model, and the development of a proper testing phase of the hazard model to quantify the consistency and skill of the hazard model. We review the state-of-the-art of the national seismic hazard modeling, and how the most innovative approaches try to address future challenges.

The Philippine archipalego is tectonically complex and seismically hazardous, yet few seismic hazard assessments have provided national coverage. This paper presents an updated probabilistic seismic hazard analysis for the nation. Active shallow crustal seismicity is modeled by faults and gridded point sources accounting for spatially variable occurrence rates. Subduction interfaces are modelled with faults of complex geometry. Intraslab seismicity is modeled by ruptures filling the slab volume. Source geometries and earthquake rates are derived from seismicity catalogs, geophysical datasets, and historic-to-paleoseismic constraints on fault slip rates. The ground motion characterization includes models designed for global use, with partial constraint by residual analysis. Shallow crustal faulting near metropolitan Manila, Davao, and Cebu dominates shaking hazard. In a few places, peak ground acceleration with 10% probability of exceedance in 50 years on rock reaches 1.0 g. The results of this study may assist in calculating the design base shear in the National Structural Code of the Philippines.

<p>As part of the 2018 National Seismic Hazard Assessment (NSHA), we compiled the geographic information system (GIS) dataset to enable end-users to view and interrogate the NSHA18 outputs on a spatially enabled platform. It is intended to ensure the NSHA18 outputs are openly available, discoverable and accessible to both internal and external users. <p>This geospatial product is derived from the dataset generated through the development of the NSHA18 and contains uniform probability hazard maps for a 10% and 2% chance of exceedance in 50 years. These maps are calculated for peak ground acceleration (PGA) and a range of response spectral periods, Sa(T), for T = 0.1, 0.2, 0.3, 0.5, 1.0, 2.0 and 4.0 s. Additionally, hazard curves for each ground-motion intensity measure as well as uniform hazard spectra at the nominated exceedance probabilities are calculated for key localities.